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SUMMARY 


This  paper  gives  a  brief  consideration  of  rotary  wing  finite  difference 
methods.  The  main  concern  is  the  specification  of  boundary  conditions  to  prop¬ 
erly  account  for  the  effect  of  the  wake  on  the  blade.  Examples  are  given  of  an 
approach  where  wake  effects  are  introduced  by  specifying  an  equivalent  angle  of 
attack.  An  alternate  approach  is  also  given  where  discrete  vortices  are  intro¬ 
duced  into  the  finite  difference  grid.  The  resulting  computations  of  hovering 
and  high  advance-ratio  cases  compare  well  with  experiment.  Some  consideration 
is  also  given  to  the  modeling  of  low-to-moderate  advance-ratio  flows.^ 
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1.  Introduction 


The.  aerodynamic  prediction  tools  available  for  routine  rotor  flow  predic¬ 
tion  and  design  have  been  limited  to  the  simplest  methods  by  the  cost  of  compu¬ 
tation.  However,  current  rapid  decreases  in  the  cost  of  computation  are  enabling 
a  quiet,  but  fundamental  reappraisal  of  our  flow-prediction  methodology.  At 
present.  Integral  methods  (i.e.,  panel  or  lifting  surface)  with  prescribed  wakes 
are  the  most  comprehensive  methods  being  used  in  a  routine  fashion.  However, 
these  methods  are  currently  unable  to  treat  nonlinear  flow  fields  —  most  notably 
those  associated  with  transonic  flows.  The  best  current  means  of  treating  such 
problems  is  with  finite  difference  methods  (FDM). 

The  essence  of  the  FDM  approach  is  to  use  Taylor  series  expansions  to 
approximate  the  partial  differential  equations  by  difference  equations,  which  are 
then  solved  by  implementing  a  variety  of  matrix  methods.  The  method  solves  the 
flow  region  at  a  finite  number  of  grid  points.  Proper  flow  resolution  depends 
on  placing  a  sufficient  number  of  grid  points  at  all  significant  flow  regions. 
Unfortunately,  a  typical  real  rotor  flow  often  contains  so  many  important  flow 
regions  (the  blade,  the  wake  system,  and  the  fuselage  and  acoustic  propagation 
regions)  that  totally  self-contained  FDM  solutions  of  complete  rotor  flows  is 
not  likely  within  the  next  five  years.  However,  practical  rotor  loads,  perfor¬ 
mance,  and  design  work  using  FDM  may  be  accomplished  using  isolated  blade  compu¬ 
tations,  if  we  know  how  to  properly  specify  our  boundary  conditions.  The  intent 
of  this  paper  is  to  demonstrate  various  types  of  boundary  condition  specifications 
to  perform  FDM  computations  of  various  rotor  flow  conditions.  For  convenience,  a 
small-disturbance  potential  equation  will  be  used  for  this  modeling. 


Equation  of  Motion  and  Boundary  Conditions 


The  simplest  model  equation  which  encompasses  transonic  flow  effects  Is 
the  classical  transonic  small-disturbance  equation,  which  for  a  rotor  has  the 
form. 


A*tt  +  B*xt 


t>  +  C(J»  +  D<f> 
rzz  yy  xy 


where  F  is  a  nonlinear  flux  term  and  is  a  function  of  <j>x  and  4>t.  This  equa¬ 
tion  was  derived  under  the  assumption  of  small  thickness  and  crossflow  and  high 
Mach  number  [1],  By  the  means  of  various  mixed  difference  schemes,  Eq.  (1)  has 
been  shown  to  give  good  transonic  solutions  for  unswept,  nonlifting  rotors  for 
advance  ratios  as  high  as  0.6  [2,3]. 


Equation  (1)  can  be  solved  in  three  forms.  1)  The  steady  form  in  which 
all  time-derivative  terms  are  ignored  [4].  In  this  case,  high-speed  hover  solu¬ 
tions  have  been  obtained  using  a  classical  relaxation  scheme  to  solve  the  differ¬ 
ence  equations.  In  ref.  [4]  the  effect  of  the  wake  was  obtained  by  including 
one  vortex  in  the  solution  grid.  This  gave  a  solution  which  bore  qualitative 
resemblance  to  a  hover  flow,  but  the  model  was  too  simple  to  be  very  accurate. 

2)  The  low-frequency  solution  in  which  all  unsteady  terms  except  4>xt  are 
ignored  [3].  This  approximation  is  valid  for  flow  variations  which  occur  in 
several  (say,  4)  or  more  chord  lengths  of  blade  travel,  which  includes  most 
rotor  phenomena.  Nonlifting  computations  using  this  model  have  been  made  using 
an  alternating  direction  implicit  (ADI)  scheme  and  have  compared  very  well  with 
experimental  data.  However,  good  lifting  solutions  have  not  yet  been  obtained 
due  to  the  lack  of  inclusion  of  a  downwash  model.  3)  The  fully  unsteady  case 
in  which  all  unsteady  terms  are  included.  This  case  includes  all  high-frequency 
flows,  the  most  notable  of  which  is  that  resulting  from  the  near  passage  of  a 


vortex.  Perhaps  the  greatest  hindrance  to  this  problem  is  the  lack  of  suitable 
experimental  data. 


In  general,  Eq.  (1)  can  be  solved  for  simple  flows;  that  is,  the  physics 
of  compressibility  no  longer  is  a  problem.  However,  the  physical  models  required 
to  accurately  simulate  real  rotor  flows  (by  FDM)  have  not  been  completely  worked 
out  or  used.  The  rotor  operates  in  the  vicinity  of  a  complex  and  often  physi¬ 
cally  close  vortex  system.  Because  we  can  only  practically  solve  for  the  flow 
on  an  Isolated  blade,  the  vortex  system  (as  well  as  other  flow  features)  must  be 
accounted  for  by  a  careful  specification  of  the  boundary  conditions.  The 
simplest  possible  set  of  boundary  conditions  for  an  isolated  lifting  rotor  must 
contain  the  following  features  (Fig.  I): 

1.  The  surface  tangency  condition.  For  the  small  disturbance  equation 
this  consists  of  equating  the  upper  and  lower  surface  slopes  of  the  blade  to  the 
vertical  flow  component  at  a  mean  surface. 

2.  The  trailing  vorticlty  sheet.  This  is  an  interior  flow  boundary 
condition  —  a  sheet  extending  rearward  from  the  blade,  across  which  there  is  a 
specified  potential  discontinuity  whose  strength  is  the  vorticity  that  has  been 
convected  from  the  trailing  edge.  This  sheet  is  usually  assumed  to  be  undis¬ 
torted  and  coincident  with  a  coordinate  surface. 

3.  The  far-field  boundary  condition.  The  outer  boundary  surfaces  are 
usually  located  sufficiently  far  from  the  rotor  that  free-stream  conditions 

(j>  =  0  or  d>n  =  0  can  be  specified  there.  It  is  important  that  the  boundary 
condition  used  on  an  outer  surface  be  compatible  with  any  interior  boundary  that 
extends  to  this  surface.  For  instance,  compatibility  with  the  trailing  vorticity 
sheet  requires  that  free-stream  pressure  be  specified  on  the  rear  boundary. 

For  a  steady  flow  this  corresponds  to  setting  <J>n  =  0. 

4.  The  inner  radial  surface  boundary.  To  conserve  mesh  points,  only  the 
outer  Dortion  of  the  rotor  flow  is  computed.  Since  the  inner  rotor  flow  is  two 
dimensional  away  from  the  tip,  it  suffices  to  use  strip  theory  (i.e.,  4>nn  =  0) 
on  the  inboard  computational  plane. 

Although  the  above  boundary  conditions  are  all  that  are  required  for  fixed 
wing  work,  they  are  incomplete  for  the  rotor  situation  because  the  downwash 
induced  by  the  wake  system  is  not  accounted  for  (except  for  that  flow  induced  by 
the  vortex  sheet  included  in  the  grid).  The  simplest  possible  way  to  include 
this  downwash  is  to  modify  the  surface  tangency  condition  by  prescribing  an 
effective  angle  of  attack.  An  example  of  this  approach  is  given  in  the  follow¬ 
ing  section. 

3.  Computation  of  High-Speed  Advancing  Blade  Loads.  Comparison  with  Experiment. 

An  ideal  situation  in  which  to  employ  the  angle-of-attack  approach  is  for 
high  advance-ratio  flight  (say.  greater  than  0.3)  where  the  induced  downweih  is 
small.  Furthermore,  there  exists  a  good  body  of  model  blade  surface  pressure 
data  for  this  case  [2,3,5].  This  data  (obtained  at  ONERA  using  2  and  3  bladed 
rotors  with  Instrumented  removable  tips)  covers  a  number  of  high  advance-ratio 
conditions  ranging  from  nonlifting  to  moderately  loaded  (Cp/cr  -  0.06)  cases. 

All  previous  computations  of  these  cases  have  been  only  for  the  nonlifting  c.^ses 
using  the  low-frequency  form  of  Eq.  (1).  To  perform  some  first  computations  of 
the  lifting  cases,  the  induced  incidence  has  been  calculated  using  the  Drees 
downwash  model  [6].  This  simple  model  neglects  all  discrete  blade  effects  (such 


as  the  tip  vortices)  and  computes  downwash  only  at  the  rotor  center  and  3/4 
radius  (the  downwash  elsewhere  is  taken  as  a  Ixnear  interpolation  or  extrapola¬ 
tion  from  these  two  points).  The  simplifications  of  this  model  are  such  that  it 
was  not  deemed  necessary  to  modify  this  downwash  by  the  vortex-sheet-induced 
velocities  when  it  was  used  in  the  finite  difference  code. 

The  unsteady  calculations  have  been  performed  on  an  isolated  blade  with 
an  angle  of  attack  prescribed  along  the  blade  for  each  azimuthal  location  by: 

a(r/R,40  -  0C  +  0W (r /R  -  0.7)  +  ax  +  a2  +  a3 


with 


ax  =  tan-1[V  sin  aq  cos  $/(0r  +  V  cos  aq  sin  ip)] 

a2  *  tan-1{  (-r0(dg/dij>)  -  V  cos  cos  \j)  sin  6] /(Or  +  V  cos  aq  sin  ip)  ] } 
a3  =  tan^f-V^  cos  3/(ftr  +  V  cos  sin  ip)  ] 
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The  experimental  values  are  used  for  the  rotor  shaft  only  (Oq),  the  col¬ 
lective  pitch  angle  (0C),  the  blade  twist  angle  (6W),  and  the  flapping  angle  (3). 
An  example  showing  the  evolution  and  relative  magnitude  of  the  different  terms: 

0C  +  0w(r/R  -  0.7);  Oj ;  a2;  and  a}  is  given  in  Fig.  2.  It  is  seen  in  this 
figure  that  although  the  Induced  angle  is  the  smallest  single  contributer  to  the 
inflow,  it  is  nevertheless  significant  in  the  first  quadrant  where  the  other 
terms  nearly  cancel  each  other. 

Two  experimental  cases  have  been  computed  corresponding  to  the  following 
conditions: 
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First  case: 


CT/o  -  0.05,  V  «  76  m/sec,  OR  *  210  m/s 


Second  case: 


C^,/o  ■  0.0665,  V  ■  81.4  m/s,  OR  ■  210  m/s 

The  theoretical  pressure  distributions  for  the  azimuths  ip  =  30°,  60°,  90°, 
120°,  150°,  and  180°  are  compared  with  the  experimental  results  in  Figs.  3-5 
(Gp/o  ■  0.05)  and  Figs.  6-8  (Cp/o  =  0.0665).  The  evolution  of  the  pressure 
distribution  is  fairly  well  predicted,  including  the  development  of  a  shock  wave 
for  ip  =  90°  and  120°.  For  the  second  case,  the  computed  shock  wave  tends  to  be 
stronger  and  located  further  aft  of  the  experimentally  obtained  shock. 

In  general,  the  discrepancies  between  experiment  and  calculation  tend  to 
be  greatest  in  the  first  quadrant.  These  differences  between  experiment  and 
calculation  are  more  easily  seen  if  we  display  the  time  history  of  the  pressure 
coefficient  at  several  points  on  the  blade  surface.  Such  time  histories  are 
shown  for  the  spanwise  section,  r/R  =  0.9,  for  Cp/o  =  0.5,  and  for  0.665 
(Figs.  9  and  10,  respectively).  The  differences  seen  in  the  first  quadrant  are 
almost  certainly  a  result  of  inaccuracies  in  the  very  simple  inflow  model. 

(Recall  that  the  induced  velocity  is  much  greater  in  the  first  quadrant  than  in 
the  second.)  The  most  significant  error  in  the  second  quadrant  is  in  the  shock 
location.  This  may  also  be  a  residual  effect  of  the  induced  velocity.  However, 
an  equally  likely  possibility  is  that  we  are  seeing  the  effects  of  neglecting 
the  boundary  layer.  It  has  been  shown,  by  a  number  of  calculations  with  boundary 
layer  models  included,  that  the  effect  of  viscosity  is  to  weaken  the  shock  and 
move  it  upstream  [7]. 

Overall,  the  results  obtained  by  these  first  lifting  calculations  of  the 
flow  around  a  helicopter  blade  in  forward  flight  (p  >  0.3)  are  in  good  agreement 
with  the  experimental  results.  However,  the  use  of  a  more  sophisticated  model 
for  the  determination  of  the  induced  velocity  is  certainly  necessary  for  the 
tower  advance  ratios. 

4.  Lifting-Rotor  Computations  by  Intra-Grid  Vortex  Insertion.  Steady  Cases. 

Transonic  computations  using  an  angle-of-attack  condition  cannot  be  valid 
when  the  vortex  is  too  close  to  the  blade  surface.  This  is  because  the  flow 
nonlinearity  demands  that  the  effect  of  the  vortex  be  computed  at  all  points  in 
the  flow,  and  not  just  at  the  blade  surface.  To  handle  this  situation,  one  can 
conceive  of  solving  for  the  perturbation  potential  about  the  vortex-induced 
velocity  field.  However,  a  simpler  approach  is  to  modify  the  existing  vortex 
sheet  logic  referred  to  in  Section  2.  That  is,  a  vortex  line  (of  constant 
strength  for  the  hover  case)  can  be  inserted  into  the  grid  by  specifying  a  sur¬ 
face  across  which  there  is  a  potential  discontinuity  whose  magnitude  equals  the 
desired  vortex  strength.  The  edge  of  this  sheet  is  the  desired  vortex  line. 

For  a  constant-strength  vortex  the  sheet  is  merely  a  branch  cut  and  its  orienta¬ 
tion  should  have  no  effect  on  the  solution.  Such  a  vortex  is  shown  in  Fig.  1, 
represented  by  a  vertical  sheet. 

Hover  flight  is  a  good  case  to  test  the  approach  described  previously 
because  the  vortex  strength  is  constant.  Furthermore,  the  vortex  geometry  is 
fairly  simple  and  can  be  reasonably  represented  by  straight  lines  parallel  to 
the  flow  in  the  near-blade  region.  For  present  purposes,  a  number  of  steady 


finite  difference  computations  have  been  made  of  hover  flows  using  the  relaxation 
code  of  Ref.  [4].  In  the  current  computations,  two  vortices  (represented  by 
vertically  oriented  sheets)  are  located  beneath  the  blade.  However,  this  does 
not  yet  account  for  the  total  effect  of  the  entire  wake  system.  We  assumed 
therefore  that  the  more  distant  elements  of  the  wake  system  can  be  accounted  for 
by  a  modified  angle  of  attack.  The  modification  consists  of  computing  the  angle 
of  attack  (using  the  Biot-Savart  law)  only  from  those  elements  of  the  wake  system 
not  contained  in  the  finite  difference  grid.  Figure  11  shows  the  computed  load 
distributions  on  hypothetical  twisted  and  untwisted  methods.  The  FDM  results  are 
compared  with  those  of  the  lifting-surface  code,  HOVER,  described  in  Ref.  [8]. 

Both  codes  use  the  identical  far-wake  model  [8]  and  assume  the  vortex  strength  to 
equal  the  maximum  blade-bound  vorticity.  The  comparison  of  the  two  methods  is 
quite  close.  The  differences  that  do  exist  could  stem  from  a  number  of  causes 
including  vortex  curvature  and  boundary  conditions. 

There  also  exists  a  convenient  body  of  experimental  data  [9]  by  which  to 
test  the  previously  mentioned  wake  treatment.  This  data  consists  of  extensive 
surface  pressure  distributions  and  vortex  geometry  measurements  for  a  high-speed 
hovering  rotor  (untwisted  blade,  NACA  0012,  /R=6.0).  The  vortex  and  blade 
geometry  is  the  only  required  input  for  computing  the  rotor  flow.  A  comparison 
of  the  experimental  lift  distribution  with  the  finite  difference  and  lifting- 
surface  methods  is  shown  in  Fig.  12.  At  this  rotor  speed  (M^  =  0.6),  the 
lifting-surface  and  finite  difference  results  compare  closely.  However,  there  is 
a  greater  difference  near  the  tip  than  occurs  in  the  previous  high  aspect-ratio 
results.  This  is  almost  certainly  a  result  of  there  being  no  wake  curvature  in 
the  finite  difference  grid.  Both  the  finite  difference  and  lifting-surface 
results  overestimate  the  maximum  lift.  Overall,  however,  the  comparison  with  the 
data  is  good.  Also  included  in  Fig.  12  is  a  finite  difference  computation  using 
only  an  angle-of-attack  boundary  condition.  For  this  computation  the  angle  of 
attack  was  obtained  from  momentum  theory.  It  is  not  surprising  to  find  that  the 
resulting  lift  is  seriously  in  error  at  the  tip,  due  to  the  lack  of  discrete  blade 
effects  in  the  momentum  model.  However,  the  inboard  results  are  good  for  this 
case.  This  indicates  that  such  a  simple  downwash  model  may  be  generally  useful 
for  specifying  the  inboard  flow  regions. 

The  comparison  of  lifting  surface,  FDM  results,  and  experiment  is  espe¬ 
cially  interesting  at  high  tip  speeds.  Figure  13  shows  such  a  comparison  at  a 
tip  speed  of  M-j  =  0.877.  Since  the  lifting-surface  method  is  linear  (as  are 
almost  all  integral  methods)  it  is  not  surprising  that  this  method  greatly  over¬ 
estimates  the  lift  in  the  tip  region.  However,  the  FDM  computations  compare 
well  with  the  experimental  lift  distribution.  Figure  13  also  shows  a  comparison 
of  the  measured  pressure  distributions  with  FDM  results.  The  comparisons  are 
favorable  even  to  the  extent  of  properly  predicting  the  correct  shock  locations. 
This  is  a  little  surprising  when  one  considers  that  the  computed  pressure  per¬ 
turbations  exceed  the  assumed  limits  for  small  perturbation  theory.  Indeed,  the 
inaccuracies  seen  in  the  leading-edge  region  undoubtedly  stem  from  small  distur¬ 
bance  errors.  In  addition,  no  boundary  layer  corrections  have  been  made. 

All  operational  helicopters  operate  at  low  enough  tip  speeds  that  linear 
integral  methods  will  always  suffice  for  hover  work.  However,  in  demonstrating 
the  ability  to  easily  duplicate  integral  results  by  FDM  and  to  effectively  apply 
these  to  high  speed,  we  will  have  performed  an  essential  validation  in  prepara¬ 
tion  for  general  forward-flight  computations. 


5.  Vortex-Induced  Effect  in  Forward  Flight. 


The  extension  of  the  previous  mode  of  vortex  representation  to  unsteady 
cases  has  geometrical,  numerical,  and  practical  problems  associated  with  it.  To 
begin,  the  wake  for  these  cases  can  no  longer  be  represented  by  a  straight, 
unmoving  vortex  of  constant  strength.  The  sheet  defining  the  vortex  system  is 
now  a  physically  real  vorticity  sheet  (for  a  tip  vortex  of  varying  strength)  and 
cannot  be  given  the  arbitrary  treatment  of  the  hover  cases.  The  tip  vortex 
assumes  a  wide  range  of  orientations  with  respect  to  the  blade  and  can  be  nearly 
parallel  to  it  in  certain  encounters.  Clearly,  a  general  forward-flight  flow 
treatment  requires  the  fully  unsteady  form  of  Eq.  (1).  Although  there  are 
methods  to  solve  the  fully  unsteady  Eq.  (1),  there  have  not  yet  been  many  experi¬ 
ments  performed  to  verify  these  methods.  A  considerable  problem  is  that  there 
currently  exists  little  usable  wake  and  loading  data  by  which  to  evaluate  new 
forward-flight  methods  at  the  lower  advance  ratios  (say,  below  0.3).  Therefore, 
at  present  the  best  approach  for  computation  work  is  to  calculate  model  problems 
which  contain  the  qualitative  effects  of  blade/vortex  interactions  in  the  simplest 
possible  way. 

One  such  model  is  the  two-dimensional  interaction  of  an  airfoil  with  a 
moving  vortex.  This  model  is  particularly  easy  to  implement  because  the  vortex 
can  (as  in  the  hover  case)  be  represented  by  an  arbitrarily  oriented  sheet 
across  which  a  potential  discontinuity  is  specified  (Fig.  14).  The  edge  of  the 
sheet  (the  vortex)  is  currently  stepped  through  the  grid  in  a  straight  line  at 
the  speed  of  the  undisturbed  free-stream  velocity.  The  strength  of  the  vortex 
is  given  as  an  effective  lift  coefficient,  CLy.  of  an  airfoil  having  the  same 
circulation  as  the  vortex.  The  technique  for  solving  the  fully  unsteady  Eq.  (1) 
in  its  two-dimensional  form  is  a  simple  ADI  scheme  given  in  Ref.  [10].  The 
effect  of  unsteadiness  in  these  computations  is  shown  in  Fig.  15,  which  compares 
blade  surface  pressure  distributions  for  a  fixed  and  moving  vortex.  These  dis¬ 
tributions  are  shown  when  the  vortex  is  at  the  leading-edge,  mid-chord,  and 
trailing-edge  chordwise  locations  and  at  0.96  chords  below  the  blade  (CLV  =  0.1). 
Although  the  steady  case  clearly  shows  the  effect  of  the  vorticity  there  is  no 
apparent  disturbance  for  the  unsteady  case.  This  result  undoubtedly  reflects 
the  fact  that  the  vortex  exerts  no  net  force  on  the  fluid  in  the  unsteady  case. 

In  fact,  in  order  to  get  a  sizable  unsteady  effect  on  the  surface  pressures  when 
the  vortex  is  beneath  the  profile,  it  is  necessary  to  bring  the  vortex  much 
closer  to  the  blade  and  increase  its  strength.  Figure  16  shows  the  lift  varia¬ 
tion  for  a  blade  with  a  vortex  of  strength  Clv  *  0.4  whose  path  is  0.26  chords 
beneath  the  blade.  As  the  vortex  passes  beneath  the  blade  it  is  seen  to  signifi¬ 
cantly  enlarge  and  strengthen  the  supercritical  flow  region  on  the  bottom  surface. 
However,  the  effect  quickly  disappears  after  the  vortex  has  passed  the  shock.  It 
should  be  recognized  that  this  latter  computation  could  contain  a  serious  error 
because  it  does  not  currently  allow  the  vortex  to  follow  the  flow  (which  varies 
greatly  from  the  free  stream  close  to  the  airfoil).  This  error  should  cause  an 
exaggeration  of  the  effect  of  the  vortex  on  the  airfoil  loads.  It  seems  likely 
therefore  that  these  close  high-frequency  interactions  will  have  a  far  smaller 
effect  on  the  loads  than  quasi-steady  theory  would  indicate.  This  does  not  say, 
however,  that  the  effect  will  be  unimportant.  The  determination  of  the  magnitude 
of  this  blade-vortex  interaction  must  await  a  good  experiment  and  more  detailed 
theoretical  modeling. 

5.  Concluding  Remarks 

The  increasing  availability  of  good  FDM  rotor  codes  and  cheaper  computers 
will  ultimately  enable  the  routine  simulation  of  the  nonlinear  unsteady 


three-dimensional  flow  about  rotors.  Initially,  these  simulations  will  be  iso¬ 
lated  blade  computations  whose  specified  boundary  conditions  account  for  the 
effect  of  the  wake  and  interference  effects.  This  paper  demonstrates  two 
approaches  to  the  treatment  of  boundary  conditions  which  result  in  realistic  FDM 
computations  of  both  hover  and  high  advance-ratio  (p  >  0.3)  rotor  flows.  For  the 
latter  case  (computed  with  a  low-frequency  unsteady  code)  discrete  vortex 
effects  are  specified  as  an  angle-of-attack  variation.  The  inflow  model  is  very 
simple  and  can  be  improved,  although  good  comparisons  with  data  have  been 
obtained.  Even  better  comparisons  have  been  made  for  high-speed  hover  computa¬ 
tions  (using  a  steady  code).  However,  for  this  case  it  is  necessary  to  use  a 
more  exact  boundary  condition  by  actually  inserting  some  of  the  tip  vortices 
into  the  finite  difference  grid.  In  addition,  a  modified  angle  of  attack  (com¬ 
puted  from  those  wake  elements  not  included  in  the  grid)  is  also  required.  For 
the  low-to-intermediate  range  of  advance  ratios,  the  problem  is  more  difficult 
because  the  wake  has  a  complicated  structure  which  cannot  be  ignored.  Neverthe¬ 
less,  demonstration  computations  (using  a  fully  unsteady  code)  have  been  made  on 
a  two-dimensional  blade/vortex  interaction  problem  which  contains  the  basic 
numerical  (if  not  all  the  geometrical)  features  of  the  problem. 

The  numerical  techniques  used  in  this  paper  are  not  new  but  they  have  not 
yet  seen  the  full  realization  of  their  capabilities  for  rotor  loads,  performance 
prediction,  and  design  work.  The  hover  and  high  simulation",  performed  in  this 
paper  may  be  the  first  truly  realistic  FDM  rotor  simulatin  There  seems  to  be 

little  technical  reason  for  reticence  in  the  use  of  FDM  te  liques  and  we  trust 
that  the  economic  reasons  will  rapidly  diminish. 
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Fig.  3.  Comparison  of  measured  and  computed  chordwise  pressure  distribution  at 
different  azimuthal  angles,  C* j/o  ■  0.05,  r/R  ■  0.85. 


Fig.  4.  Comparison  of  measured  and  computed  chordwise  pressure  distribution  at 
different  azimuthal  angles,  G p/a  *  0.05,  r/R  -  0.90. 
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Fig.  5.  Comparison  of  measured  and  computed  chordvise  pressure  distribution  at 
different  azimuthal  angles,  Cj/o  -  0.05,  r/R  -  0.95. 
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Fig.  11.  Comparison  of  computed  loadings  by  panel  and  finite  difference  method 
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.  Comparison  of  measured  and  computed  chordwise  pressure  distribution  at 
877,  r/R  -  0.89,  and  r/R  -  0.96,  respectively. 
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Fig.  13.  Comparison  of  measured  and  computed  chordwise  pressure  distribution  at 
Hj.  ■  0.877,  r/R  -  0.5,  and  r/R  -  0.68,  respectively. 
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SMALL  DISTURBANCE  BOUNDARY  CONDITIONS 
FOR  2-D  WING/VORTEX  PROBLEM 


Fig.  14.  Boundary  conditions  for  the  blade/vortex  interaction  problem 
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Pig.  16.  Unsteady  flow  computations  for  the  blade/vortex  interaction  problem 
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